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PASCO; STRUCTURAL PANEL ANALYSIS AND SIZING 


CODE , CAPABILITY AND ANALYTICAL FOUNDATIONS 
W. Jefferson Stroud and Melvin S. Anderson 

SUMMARY 

A computer code denoted PASCO which can be used for analyz- 
ing and sizing uniaxially-stif f ened composite panels is described. 
Buckling and vibration analyses are carried out with a linked- 
plate analysis computer code denoted VIPASA, which is incorporated 
in PASCO. Sizing is based on nonlinear mathemctical programming 
techniques and employs a computer code denoted CONMIN, also 
incorporated in PASCO. Design requirements considered are initial 
buckling, material strength, stiffness, and vibration frequency. 
The report describes the capability of the PASCO computer code 
and the approach used in the structural analysis and sizing. 

INTRODUCTION 

Stiffened panels made of metal and/or composite materials 
have wide application in aerospace structures. These panels are 
generally designed to have low mass and must meet niomerous design 
requirements involving, for example, buckling, stiffness, 
material strength, and limitations on panel geometry. In an 
effort to increase the structural efficiency of these panels, 
design concepts are being explored which exhibit complex 
buckling modes, requiring sophisticated stability analyses. In 



addition, composite panels may require relatively sophisticated 
stress analyses. 

To address these needs, a computer code denoted PASCO has 
been developed for analyzing and sizing stiffened composite panels. 
The code attempts to balance the user requirements of generality, 
simplicity, rigor, and modest computer resources. This report 
describes the analytical foundations of PASCO to the extent that 
it would h ilp a user understand the analysis and sizing procedures, 
select appropriate options, and interpret answers. Complex 
theoretical discussions are treated in the references. The users 
manial for PASCO, reference 1, includes an explanation of 
structural modeling for PASCO, a discussion of program input and 
output, and several illustrative examples. Design studies 
carried out with PASCO are described in references 2 and 3 . 

Previous work dealing with the analysis and sizing of stiffened 
composite panels is discussed in references 4 and 5. 

The present report begins with a discussion of the program 
capability and approach. There follows a discussion of the 
stress and buckling analyses. Finally, the structural sizing 
strategy is discussed. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The 
calculations were made in U.S. Customary Units. In many cases, 
the FORTRAN name of the variables used in PASCO is included in 
the definition. 
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A 

A 




jk 


A., 

3k: 


panel planform area shown in figure 25 

for closed section stiffeners, area enclosed by 
the closed section in one period 

2 

A^^ - /^22 plate element i 

laminate inplane stiffnesses and smeared ortho- 
tropic inplane stiffness defined by equation (1) 

value of laminate inplane stiffness A., for 


plate element i 


jk 


AllL, AllU 
A33L, A33U 

ALLOW 

b 


lower and upper bounds on smeared orthotropic 
stiffnesses A^^ and A^^ 

allowables used in the material strength criteria 

plate element width 

width of plate element i 

v.’idth of one period of the stiffened panel 
(see figs. 6 and 8) 

quantities defined by equation (18) 

quantities associated with temperature, defined in 
equations (1) and (5) 

parameter used during sizing for prescribinc 
margin of safety on buckling load with haif- 
wavelength X equal to L/m, defined in 
equation (47) 

CONVl, C0NV2 convergence criteria for eigenvalue analysis 


^i 


C. 

1 

*^jT 


CLAM (A) 


D -1, 

3k 


D . 


3ki 


laminate bending stiffnesses and smeared ortho- 
tropic bending stiffnesses defined by 
equation (35) 

value of laminate bending stiffness D.. for 


plate element i 


jk 


DHL, DllU 
DVMOV 
e, ECC 


lower and upper bounds on smeared orthotropic 
bending stiffness 

parameters used to determine move limits during 
sizing, defined in equations (54) and (55) 

bow at panel midlength 
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E^, El 

E^, E2 

E^2' E12 

f, FREQ 
F, FACTOR 

F(y) 

Fi, F2, Fii, 
^ 12 ' ^ 22 ' ^66 

G 

GRANGE 

ITHERM 

L 

M 

m 

MAXL 

MINLAM 

NEIG(m) 


Young's modulus of composite material in fiber 
direction 

Young's modulus of composite material in direction 
transverse to fiber direction 

Shear modulus of composite material in material 
coordinate system 

frequency 

scale factor that relates the input or design 
load to the load that causes buckling or 
vibration, defined in equation (39) 

buckling displacement function, defined in 
equation (36) 

real and imaginary parts of F(y), defined in 
equation (37) 

functions that appear in Tsai-W . material 
strength criterion (see eq. (49)) 

behavioral constraint (see eq. (46) ) 

constraint deletion parameter 

parameter used to indicate the manner in which a 
bending moment produced by temperature or 
transverse load is to be treated 

panel length 

bending moment, per unit width 

applied bending moment per unit width (see figs. 

2 and 7) 

half-wavelength number, L/X 

maximum number of values of X for which buckling 
or frequency constraints are calculated 

integer that specifies smallest value of X 
for which buckling loads are examined 
(X = L/MINLAM) 

number of eigenvalues determined at X = L/m 
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Njj, Ny, Nxy, inplane longitudinal, transverse, and shear loads 
NX, NY, NXY per unit width, applied to panel 

N value of N that corresponds to eigenvalue 

^cr ^ 


N 


X 


input 


input value of N^; also value of N for 
which panel is designed 



Euler buckling load of panel 



XYi 



inplane longitudinal, transverse, and shear loads 
per unit width, applied to plate element i 



P, PRESS 


Q' 


jk 



RHO 

S 

^allow 

^i 

S 

P 

SHEAR 


t 


shear load applied to substructure 


uniform lateral pressure 

lamina stiffnesses in material coordinate system, 
defined in equations (6) through (9) 

lamina properties associated with thermal expan- 
sion; primed quantities indicate material 
coordinate system; unprimed quantities indicate 
panel coordinate system 

density 

lamina stress or strain (see eq. (48)) 

allowable value of S (see eq. (48)) 

shear flexibility for plate element i (see 
eq. (30)) 

shear flexibility for substructure p (see 
eqs. (31) and (33)) 

parameter used to indicate whether the standard 
VIPASA analysis is to be used for the X = L 
buckling load (SHEAR =0) or whether the 
adjusted analysis is to be used (SHEAR > 0) ; 
appropriate only for cases where the loading 
involves shear 

thickness 
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I 


■w 


'i 





u, V, w 


plate element displacements in local plate element 
coordinate system 

sizing variable 

X, value of X. at initial point of Taylor series 

expansion^of constraints 

X, y, z coordinate directions in local plate element 

coordinate system; axes defined in figure 6 

Zj distance from reference surface to centroid of 

plate element i 

distance from centroid of cross section to 
centroid of plate element i 

oi 2 i coefficient of thermal expansion of composite 

ALFAl, ALFA2 material in material coordinate system 

y F/F(X = L) 

Vt, TEM change in temperature 

^ v 7 laminate strains in plate element coordinate 

^ system 

6,- THET angle between material coordinate system and plate 

element coordinate system (see figure 6) 

ANUI Poisson's ratios of composite material in material 
coordinate system 

X buckling o’- vibration half-wavelength 

DESCRIPTION OF GENERAL APPROACH USED IN AND 
CAPABILITY OF PASCO COMPUTER CODE 
Overview 

PASCO has been developed to aid the engineer in the analysis 
and sizing of prismatic structures such as those shown in 
figure 1. Because of their wide application in aerospace 
structures, stiffened panels are given special emphasis in PASCO. 
For exarole, complex panel configurations can be built up from a 
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Figure 1.- Examples of typical structures, 
a relatively small number of repeating elements, the loadings 
(figure 2) available in PASCO are the type usually associated 
with panels, and practical panel desigii considerations such as 
an overall bow-type initial imperfection (figure 3) can be 
accounted for. 



Figure 2.- Available loadings Figure 3.- Overall bow-type 

applied to hat-stiffened panel. initial imperfection. 

The panel cross section is composed of an arbitrary 
assemblage of thin, flat, rectangular plate elements that are 
connected together along their longitudinal edges. Each plate 
element consists of a balanced symmetric Icuninate of any number 
of layers of orthotropic material. Any group of element widths. 
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layer thicknesses, and layer orientation angles can be selected 
as sizing variables. For example, in the blade-stiffened panel 
configuration shown in figure 4, the blade depth can be allowet 
to vary, the overall stiffener spacing can be held fixed, the 
thickness of the material at the 0° and ±45° orientations can be 
allowed to vary, and the orientation angles themselves can be 
held fixed. 



2 


Figure 4.- Blade-stiffened panel configuration. 

When used in the analysis mode, PASCO can be used to cal- 
culate laminate stiffnesses, lamina stresses and strains 
(including the effect of temperature), buckling loads, vibration 
frequencies, and overall panel stiffness. When used in the siz- 
ing mode, PASCO adjusts the sizing variables to provide a low- 
mass panel design that will carry a set of specified loadings 
without failure by buckling or material strength and that will 
meet other design requirements such as upper and lower bounds on 
sizing variables, upper and lower bounds on overall bending. 
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extensional and shear stiffness, and lower bounds on vibration 
frequencies. 


Approach and Capability 

The approach used in and the capability of PASCO are summar 
ized in figure 5. The topics listed in figure 5 are discussed 
briefly here and are explained in greater detail in subseq-"ent 
portions of the report. 


DESIGN CONDITIONS 

• N , N , N 

X y xy* X 

• Lateral pressure 

• Bow-type imperfection 

• Temperature 

• Multiple sets of design conditions 


ANALYSIS 


• VIPASA for eigenvalues 

• Prebuckling stresses include bending stresses ;aused by applied 

moment, lateral pressure, bow, temperature, and transverse load 

• Lamina stresses and strains for strength criteria 

PANEL CONFIGURATION AND MATERIAL 


f General configuration 
t Multiple materials 

• Materials orthotropic at arbitrary angle 

• Balanced, symmetric laminates 
f Uniform along panel length 

• Detailed modeling 

SIZING, SIZING VARIABLES. AND CONSTRAINTS 

• Nonlinear mathematical programming 

• Minimize mass 

§ Sizing Variables; element widths, ply thickness, and ply 
orientation angle 

• Sizing variable linking 

• Bounds on sizing variables 

• Constraints: buckling, material strength, stiffness, 

vibration frequency 

Figure 5.- PASCO capability and approach. 




Design conditionb . - ne design conditions considered are a 
loading of N , N , N , and M , lateral pressure, an over- 
all bow-type imperfection, and temperature. Lateral pressure 
and the bow are treated using a beam-column approach. Thermal 
stresses are calculated assuming that the temperature (or, more 
precisely, the change in temperature) is specified in each plate 
element. Temperature does not vary with the sizing variables. 
Panels can be sized for multiple sets of design conditions. 

Analysis . - The eigenvalue analyses for buckling and 
vibration frequency are performed by a stiffened panel analysis 
code denoted VIPASA (refs. 6 and 7) . The prebuckling stress 
state includes bending stresses caused by an applied bending 
moment, lateral pressure, an o erall bow, temperature, and a 
transverse load. Resultant prebuckling lamina stresses and 
strains are calculated for the material strength criteria. 

Panel configuration and material .- The panel cross section 
can have a general configuration. Each plate element can contain 
multiple orthotropic materials oriented at arbitrary angles. 
However, each plate element must be a balanced, symmetric 
laminate and must be uniform along its length. Many unsymmetric 
laminates can be generated by stacking plate elements composed 
of symmetric Icuninates. (This approach is discussed in refer- 
ence 1.) Curved panels can be modeled using a series of flat 
plate elements. Provision for offsets and a large ’"umber of 
distinct layers in each Iciminate allow relatively detailed 
modeling of the panel cross se''tion. 
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Sizing, sizing variables, and constraints .- Sizing is 
carried out using a nonlinear mathematical programming approach 
in which the design requirements are treated as inequality 
constraints. The objective function, the quantity that is 
minimized, is the panel mass per unit width. The panel length 
is fixed. 

Any set of plate element widths, layer thicknesses, and 
layer orientation angles can serve as the sizing variables. 

The other plate element widths, layer thicknesses, and layer 
orientation angles can be held fixed or linked linearly to those 
that serve as sizing variables. During panel sizing, linking 
can be used to provide practical proportions, calculate 
offsets that change as thicknesses change, and maintain overall 
panel width. 

The design requirements, or constraints, that can be 
specified are upper and lower bounds on sizing variables, lower 
bounds on buckling and material strength, upper and lower bounds 
on overall bending, extensional, and shear stiffnesses, and 
lower bounds on vibration frequency. Separate margins of 
safety can be placed on each buckling or vibration mode. Several 
material strength criteria are included, and, if desired, the 
user can incorporate his own material strength criterion by 
writing an additional subroutine. 
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STRUCTURAL ANALYSIS 




Analytical foundations for PASCO are discussed in terms of 
structural analysis and in terms of sizing. The focus of this 
section is structural analysis. Sizing is considered in a 
subsequent section. 

Prebuckling Stress Analysis 

The prebuckling load distribution on each plate eleiaent is 
required for the buckling analysis and is used to compute 
lamina stresses and strains for the material strength criterion. 

The loads on each plate element are calculated using the follow- 
ing approach: The load distribution is first determined for a 

uniform longitudinal strain. Additional bending loads are then 
calculated and added to the load distribution determined for 
uniform longitudinal strain. Finally, the shear stress 
distribution is computed and added. Each of these steps is 
discussed in greater detail in subsequent sections of this 
report. 

Elastic relations .- The elastic relations are presented 
for a plate element with coordinate system, displacements, and • 

loading as shown in figure 6. 

Wif* the assumption of balanced, symmetric laminates, the 
general plate constitutive equations uncouple. The equations 
for inplane loads for plate element i reduce to 
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(a) Hat-stiffened panel with 
applied loading 



(W Plate element coordinate system, 
displacements, and piy 
orientation angle 


tc) Plate element i with inolane loading. 
Directions shown are positive for 
prebuck ling loads. 


Figure 6.- Plate element coordinate system, displacements, 
loading, and sign convention. 
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in which N 


^i 
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^i' 


and N 


xy; 


are the inplane loads on plate 


element i (positive in compression); A., are laminate 
stiffnesses; 
given by 


£y/ are strains (positive in compression) 


e 


X 


3u 

3x 


( 2 ) 


e 


y 


3 v 

3 y 


( 3 ) 


Y 


xy 



( 4 ) 


u and V are prebuckling displacements; and are the 

temperature terms given by 




(VT) dz 


( 5 ) 


in which VT is the change in temperature. The change in 
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temperature VT is allowed to be ply dependent and, therefore, 
call vary wxch z. 

For an orthotropic leimina with material coordinate system 
ir ;lined at an angle 0 to the plate element coordinate system 
(.ig. 6), the are first calculated in the material 

coordinate system (the primed system) and then transformed to 
the plate element coordinate system as follows: 


Q' 
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1 - v^V2 


( 6 ) 


Q* 


^2^1 


^ 1®2 


12 


1 - v^V2 1 - v^V2 


( 7 ) 


Q’ 


22 


1 - v^V2 


( 8 ) 


Q' 


33 


E 


12 


( 9 ) 


^'iT ®'ll“l 12*^2 


Q' 2 T “ ^'l 2 °^l ^' 22^2 


Qq^iji Q cos 6 t Q 2 'j> sin 0 


( 10 ) 

( 11 ) 

( 12 ) 
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cos 


e 


(13) 


'2T 


= sin^ e + Q* 


2T 


Q-T “ ^Q'it " ®'2T^ 

Uniform longitudinal strain .- With the assumption that the 
prebuckling longitudinal strain is uniform over the panel 

cross section, the strain is given by 


e.. = 



(15) 


in which is the applied longitudinal load per unit width, 

b is the width of one period of the stiffened panel (fig. 6) , 

s 

the summation extends over all elements in a period, and 


= A^^ - (Aj^ 2^^/'^22 plate i 

b^ = width of plate i 

Ci = - C 2 ,j,)/A 22 + for plate i 

The longitudinal loading N in plate i is then given by 

^i 



(16) 

(17) 


(18) 


(19) 
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In the expression for , the transverse load in 

plate i can be determined two ways. One way is to use the 

PASCO modeling rules discussed in reference 1. The other way is 

to specify the values of N with progreim input. In general, 

^i 

the modeling rules are designed so that the full is carried 

by the skin, and no N is carried by the stiffener elements. 

y 

Note that even though the longitudinal strain is 

uniform, the term in equation (19) may produce a net 

bending moment about the centroid of the panel cross section. 

Bending loads .- In this section, expressions for bending 

strains caused by various loadings are developed. These bending 

strains are combined with the uniform strains from equation ( 15) 

to calculate the total axial loading N in each plate element. 

*i 

Bending loads can be caused by an applied bending moment, a bow- 
type imperfection, lateral pressure, temperature, and/or an 
applied transverse load. These bending loads are calculated by 
PASCO and, except for those bending loads already included in 
C^ (eg. (18)), are added to the longitudinal load distribution 
given by equation (19). Certain options involving the bending 
loads produced by temperature and/or an applied transverse load 
are discussed later in this section. 

For combinations of applied moment, lateral pressure, and 
initial bow, the maximum bending moment, which occurs at panel 
midlength, is given as in reference 8 oy 
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e 


1 


( 20 ) 


N 


M = -f 


1 - Y 


PL 


- j^sec(j/Y) - ij 


wfare 



E 


and N is the applied longitudinal load, N is the Euler 

X Xg 

buckling load for the panel, is the applied bending moment 

on the panel, e is the bow at panel midlength, P is the 
lateral pressure loading, and L is the panel length. Most of 
these quantities, together with the overall panel coordinate 
system, are shown in figure 7. Note that the bending load 
caused by the applied bending moment M is not influenced by 

A 

inplane loads . 

Because the VIPASA buckling analysis requires that the stress 
distribution be constant along the panel length, the conservative 
assumption is made that the bending moment given by equation (20) 
is the bending moment over the entire panel length. 

For buckling modes having a half-wavelength X equal to the 
panel length L, the bending moments caused by a bow and/or 
lateral pressure are omitted from equation (20) . Only the 
applied bending moment is retained. 
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K 


i 



Figure 7.- Panel with applied bending moment, initial bow, 

and lateral pressure. 

Strictly speaking, equations (20) and (21) are appropriate 
only when is the sole inplane loading; however, in PASCO, 

these equations are applied to problems with combined loads. 

For combined loads, the parameter y is defined as 

F 

Y = (22) 

F(X = L) 


in which F is a scalar defined by 




buckling or vibration 


(23) 
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The input vector on the left, which is the design loading and 
frequency requirement, is scaled up or down with the parameter 
F to obtain that combination that causes buckling or vibration.^ 
Since there is a question about the validity of equation (20) 
in the case of combined loads, and since it would be inappropriate 
to use a frequency requirement to calculate y for equation (20) , 
a user should exercise cautioii in the application of equation (20) 
to calculate bending loads. For the latter reason, subsequent 
discussions in this section will focus on buckling. 

During the buckling analysis, the buckling load is calcu- 
lated for many values of buckling half-wavelength X, and 
more than one buckling load can be calculated at a given value 
of A. There is a value of F associated with each of these 
buckling loads. In equation (22), F(A = L) is the value of 
F associated with the lowest buckling load for A = L. The 
value used for F in the numerator of equation (22) depends upon 
whether the bending moment in equation (20) is to be used for 
material strength calculations or buckling calculations. 

• For material strength calculations, the bending moment that 
is used to calculate lamina stresses and strains is based on one 
of two values of F for the numerator of equation (22) . 


more complete discussion of F 
subsequent section entitled FACTOR and 


is presented in a 

F. 


20 



(1) If F(A = L) is greater than 1.0, then F = 1.0. 

(2) If F(X = L) is equal to or less than 1.0, then F 
is the value of F for the minimum buckling load for X 
considered . 

• For buckling calculations, F that appears in the numerator 
of equation (22) is the value of F associated with the eigenvalue 
number and buckling half wavelength being examined. The resulting 
bending moment is used to calculate prebuckling plate element 
loads . 

In the discussion following equation (19) , it was pointed 
out that the term, which accounts for temperature and 

transverse loads, can produce a bending moment in the panel. 

This bending moment is treated in PASCO in one of two ways: 

(1) the panel is allowed to take on a bow (JTHERM=1) , or (2) the 
panel is forced to remain flat (ITHERM=0). If the panel is 
allowed to take on a bow, the magnitude of the bow is calculated 
to produce zero bending moment in a panel loaded only by tempera- 
ture and transverse loads. This bow is then added to any initial 
bow that exists in the panel. If the panel is forced to remain 
flat, no additional bow is added and any bending moment produced 
by the terms in equation (19) is retained. The user 

selects the desired approach with the input parameter ITHERM. 

The treatment of the bending moment caused by C. can also 
be described in the following way. Let a moment M bo defined 
by 
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(24) 


N (e + e_) PL‘ 
M = ^ 

1 - Y 




•.ec(y/Y) ■ 1| “ 


■ 


M 


in which and are the only terms that do not appear in 

equation (20) . The quantity e^ is the magnitude of a bow 

calculated to produce zero bending moment in a panel loaded only 

by temperature and transverse loads, and is the r.oment 

caused by (eq. (18)). When ITHERM is set equal to 1 with 

program input, the bending moment that is added to the stress 

state associated with uniform strain is given by equation (24) . 

When ITHERM is set equal to 0 with program input, the bending 

moment that is added to the stress state associated with uniform 

strain is given by equation (20) . 

Bending loads are applied to the panel with an loading 

^i 

that varies by steps in the z direction. An example is shown in 

figure 8. In this example, a blade-stiffened panel is subjected 

to longitudinal compression and a bending moment that puts the 

skin in additional compression. The blade is modeled as three 

separate plate elements. The moment- induced load in each 

^i 

plate element, including the skin, is calculated by assuming 
that (1) the strain at the centroid of each plate element forms 
a "linear" strain distribution, (2) the net inplane load caused 
by the bending strain is zero, and (3) the net bending moment 
produced by bending strains is equal to the calculated moment. 

The resulting bending strain distribution is given by 
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I ^r-Refereiue surface I 

I ^Element 1 / (for this example) | 



Skin 



1 

1 

Blade 

r 


_i 

1 

' 



(a) Blade-stiffened panel configuration 


I 

z 

(b) for small bending moment 


7 

(c) Ejj for large 


wment 


Figure 8.- Idealized longitudinal strain distributions 
on compression panels with bending moment. 


1 


(25) 


in which is the bending strain in element i, is the 

”i 

distance from the reference surface to the centroid of plate 
element i (fig. 8) , and 


A B 


y)A.b. Z. 
^111 

X)A.b. 

^11 


(26) 


"b ' 


EAib,(Z,) 


EVi 


(27) 
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for plate element 1 


(28) 


A. = A^^ - (A,^)^'A 


12' 


'22 


The summations are over all elements in one j-eriod of the stiff- 
ened panel, and b is the width of one period. 

The bending strains from equation (25) are combined 

i 

with the uniform longitudinal strain from equation (15) 

to produce the asultant longitudinal loading in each element 


N 


X 


i 



+ 


e 



+ 


^i 


(29) 


Shear stress .- The shear stress in each plate element is 

calculated using a generalization of the approach of reference 5 

in which equilibrium and compatibility of displacement are 

employed. Define a shear flexibility 'or a plate of width 

b. 

1 


Si = bi/Aji^ 


(30) 


For plates or substructures connected in series, (^ig. 9) , 

the shear flexibility S for the substructure is given by 

P 


S = S-, + + . . . + s 

p 1 2 n 


(31) 
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^T) PKk cttiwflt or substructure nuiter 

Figure 9.- Plate elements or substructures connected in series, 
and the shear in each plate or substructure is 


N = N 
xYi XV2 


. = N = N 


in which N is the shear load on the substructure. 

For plates or siibstructures connected in parallel (fig. 10), 
the shear flexibility is given by 


1/S^ * 1/S, + 1/S, + . . . + 1/S 
p X z n 




© 

©<: / 




Plate element or substructure number 


Figure 10.- Plate elements or substructures connected in parallel. 



» 


and 





( 34 ) 


In PASCO, the overall panel shear stiffness is calculated 

by VIPASA. However, the N and S. given above could be 

xy^ 

used to calculate a total shear angle which, in turn, could be 
used to define an overall panel shear stiffness. Usually, the two 
approaches give results in close agreement. 

The shear stress in any element can also be specified 
by the user with program input. 

VIPASA Buckling Analysis 

Buckling and vibration analyses in PASCO are carried out 
with a stiffened panel analysis code denoted VIPASA (Vibration 
and Instability of Plate Assemblies including Shear and 
Anisotropy) described in references 6 and 7. For simplicity, 
only buckling terminology is used in the present discussion 
of VIPASA. It is understood, however, that the discussion 
also applies to the vibration analysis. 

The VIPASA analysis treats an arbitrary assemblage of plate 
elements with each plate element i loaded by N # ^ 

*i Yj^ 

. The buckling analysis connects these individual plate 
elements and maintains continuity of the buckle pattern across 
the intersection of neighboring plate elements. Several 
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buckling modes are shown in figures 11 and 12. VIPASA considers 
only initial buckling. Postbuckling response is not considered 
by VIPASA or PASCO. 

OVEfUa BUCKLING MODE 



BUCKLINC-MOOE SHAPE 
UNDEFOBMED SHAPE 


LOCAL BUCKLING MODE 

Figure 11.- Typical buckling modes for hat-stiffened panel. 




BUCKLING-MOPE SHAPE 
UNOEFORMEO SHAPE 


Figure 12.- Typical buckling mode for corrugated panel. 

Elastic relations .- In the VIPASA analysis, a local 
coordinate system is defined for each individual plate element. 

In the example shown in figure 6, the X, Y, Z axes define 
the local coordinate system in the local longitudinal, transverse, 
and lateral directions, respectively. The buckling displacements 
u, V, and w are defined in this local coordinate system and 
are the Scune as those shown in figure 6. 

During buckling, the out-of-plane elastic deformations of 
each plate element are defined by 
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(35) 


where M , M , and M are perturbation bending and 
^i ^i *^i 

twisting moments per unit length on plate element i, D.. are 

the Iciminate stiffnesses, and w is the perturbation displacement 

in the z-direction. Because of the requirement that each plate 

2 

element consist of a balanced, symmetric laminate, anisotropic 
effects are limited to those produced by the ^23 

terms. Therefore, in subsequent discussions, anisotropy refers 
only to anisotropy in the bending stiffness. 


5 

VIPASA does not require that plate element laminates be 
balanced and symmetric. However, as can be seen from equations (1) 
and (35), VIPASA ignores extension- she ar coupling and membrane- 
bending coupling in each plate element. The resulting elastic 
relations are the Scune as those that are obtained for a laminate 
that is balanced and symmetric. Because of the elastic relations 
in VIPASA and because balanced and symmetric laminates are the 
most common laminates in aerospace applications, PASCO input 
provides only balanced, symmetric laminates for each plate element. 
By stacking symmetric laminates (ref. 1) , many unsymmetric 
laminates can be modeled and the coupling action in the 
elastic response for these unsymmetric laminates can be 
accounted for. 
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Buckling displacements and boundary conditions .- The 
buckling displacement w assumed in VIPASA for each plate 
element is 

w Re ^F(y) (36) 

with similar expressions assumed for the inplanc* displacemei.^^ 
u and V. For F(y) written as 

F(y) = fj^(y) + i^2^y^ 

the buckling displacement w can be written as 

w = fj^(y) cos ^ - f 2 (y) sin ^ (38) 

Neglecting boundary conditions, the displacement shape assumed 
in equation (38) provides an exact solution to the governing 
differential equations if the panel and loading are uniform 
in the x-direction. The governing equations are based on the 
Kirchof f-Love hypothesis applied to each plate element. 

The functions f^^(y) and ^2^^^ allow various boundary 
conditions to be prescribed on the lateral edges of the panel. 
These boundary conditions, which include free, simple support, 
clamped, and symmetry, are discussed in the users manual, 
reference 1. Boundary conditions cannot, however, be prescribed 
on the ends of the panel. 
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Orthotropic panels with no shear loading .- For orthotropic 
panels with no shear loading, £2# the imaginary part of 
F(y), is zero. The solution fj^(y) cos provides a series 
of node lines that are straight, perpendicular to the longitudinal 
panel axis, and spaced X apart as shown in figure 13. Along 
each of these node lines, the buckling displacements satisfy 
the following simple support boundary conditions: u is 

unrestrained, v = w = 0, and w „ is unrestrained. For 

9 X 

values of X given by X = L, L/2, L/3, . . . L/m, where L 
is the panel length and m is an integer, the nodal pattern 
shown in figure 13 provides simple support boundary conditions 
at the ends of a finite, rectangular panel. An exeunple in 
which X = L/2 is shown in figure 14. 
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Figure 13.- Node lines produced by w = fi(y) cos -j- 
for orthotropic panels with no shear loading. 
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Figure 14.- Buckling of orthotropic panel under longitudinal 
loading. Mode shovm is m = 2. 

Anisotropic panels and/or panels with a shear loading .- For 

anisotropic panels and/or panels with a shear loading, £ 2 ^ 0 , 

The functions f and f 2 are such that node lines are skewed 

and not straight, but the node lines are still spaced A apart 

as shown in figure 15. In this case, the solution given by 

equation (38) is accurate only when many buckles form along the 

panel length, in which case boundary conditions at the ends are 

not important. An example in which X == L/4 is shown in 

figure 16. 

As A approaches L, the VIPASA buckling analysis for a 

panel loaded by N can be quite conservative. One explanation 

xy 

is as follows: As can be seen in figure 16, the skewed nodal 

lines given by VIPASA in the case of shear and/or anisotropy do 
not coincide with the end edges. Forcing node lines to 
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Figure 15.- Node lines produced by w = f^(y) cos ^ 


f2(y) 


for anisotropic panels and/or panels with a loading that 
includes shear. 


sin 
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Figure 16.- Buckling of panel under shear loading. 
Mode shown in m = 4. 



coincide with the en«. edges produces long-wavelength buckling 
loads that are, in many cases, appreciably higher than those 
determined by VIPASA. Calculations have shown that for long- 
wavelength buckling modes, the effect of anisotropy is minimal 
lor most practical cases. Anisotropy therefore, causes 
negligible conservatism in a VIPASA analysis. The presence of a 
shear loading can, however, lead to very conservative results 
for X equal to L. (See, for example, ref. 9.) 

Because of VIPASA 's conservatism in the case of long- 
wavelength buckling if a shear load is present, an adjusted 
shear analysis procedure can be used (at the user's option) 
for the case X = L. That adjusted analysis is discussed in a 
subsequent section entitled Adjusted Analysis for Shear 
Buckling. 

Example .- A buckling response diagram, such as that shown 
in figure 17, provides a convenient means of studying the 
buckling response of a panel and can be used to help explain 
some of the features of the buckling analysis and the computer 
code. The exeimple shown in figure 17 is for a blade-stiffened 
panel having arbitrary but reasonable proportions. The panel 
has a length L of 0.76m (30 in.) and is modeled with 
16 stiffeners. The boundary conditions on the lateral edges of 
the panel are taken to be simple support. Anisotropy is ignored. 
The loading on the panel is pure longitudinal compression; 
transverse and shear loads are taken to be zero. In the diagreun, 
the buckling load is given as a function of the 
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Figure 17.- Longitudinal buckling load as a function of 
buckling half-wavelength for blade-stiffened panel. 

nondimensional ha If -wave length X/L. The half -wavelengths 

examined by PASCO are X = L, L/2, L/3, L/4, ...., L/MINLAM 

in which MINLAM Js program input. For this excunple^ the lowest 

buckling load has a half-wavelength X = L. The buckling mode 

shape for this mode is shown in figure 18a. The next lowest 

buckling load is a relative minimum that occurs for X = L/8. 

The mode shape for this local mode is shown in figure 18b. 
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lb) LOCAL MODE . IX-U»i 


BUCKLING MODE SHAPE 

UNOEFORMEO SHAPE 

Figure 18.- Buckling mode shapes for blade-stiffened 

panel excunple. 

Although the program makes simple, exploritory calculations 
for many values of X (X = L, L/2, L/3, . . L/MINLAM) , it 
calculates the buckling load for only certain values of X. The 
program always calculates the buckling load for X = L. The 
program also calculates the buckling load for specified values of 
X given by X = L/NLAM where the vector parameter NLAM is 
input. (In order to obtain the data for figure 17, the vector 
NLAM was set equal to NLAM =2, 3, 4, 5, . . ., 30.) In addition, 
the progreun calculates any buckling load which is a relative 
minimum (X = L/8 in figure 17) that is lower than the next 
preceding calculated buckling load. Wavelengths are considered 
in order of decreasing length: L, L/2, L/3, . . . , L/MINLAM. 

Referring again to the exeunple of figure 17, if no value of NLAM 
were input, the only buckling load calculated would be for 
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A = L. The buckling load at the relative minimum X > L/8 
would not be calculated because that buckling load is greater 
than the next preceding calculated buckling load, which is at 
A = L. If, on the other hand, NLAM = 2 were input, then the 
buckling loads would be calculated for A - L, A - L/2, and 
A = L/8. The buckling load would be calculated for A « L/8 
because it is a relative minimum and because it is lower than 
the next preceding calculated buckling load - the load for 
A = L/2. 

The program input parameter NEIG(m) can be used to calculate 
more than one buckling eigenvalue at a given value of A. 

Element m in vector NEI6(m) is the number of eigenvalues 
requested at a half-wavelength of A = L/m. For example, for 
two eigenvalues at A = L, the input is NEIG(l) - 2 . In 
figure 17, the second eigenvalue at A = L is indicated by the 
square symbol at A = L. 

As explained earlier, PASCO can also account for an overall 
bow-type initial imperfection. The buckling response curves 
shown in figure 19 are for the same blade-stiffened panel discussed 
above, but with three different assumptions regarding an overall 
bow: (1) a positive bow of e/L = +0.003, (2) a negative bow of 

e/L = -0.003, and (3) a zero bow, e/L = 0.0. As in figure 17, the 
only loading is longitudinal compression. The curve for 
e/L - 0.0 is the same as that shown in figure 17. The bow does 
not directly affect the buckling load for A = L. For this 
reason, the panel has the same buckling load at A » L for the 
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Figure 19.- Longitudinal buckling load as a function of 
buckling half-wavelength for blade-stiffened panel with 
positive bow, negative bow, and zero bow. 

positive, negative, or zero bow. For a positive bow, which causes 

additional compression in the skin, the lowest buckling load 

occurs at X » l/ 30. For a negative bow, which causes additional 

compression in the tip of the blade, the lowest buckling load 

occurs for X L/8. 

FACTOR and F 

In VIPASA, FACTOR is the unknown in the eigenvalue analysis. 
The desired eigenvalue is identified by half-wavelength X and 
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by the eigenvalue number at thac value of X. A buckling analysis 
in VIPASA is merely an eigenvalue analysis at zero frequency. 

The eigenvalue solution technique in VIPASA can bo summarized 
as follows. For any set of values of FACTOR and ha If -wave length 
mathematical expressions in VIPASA provide the number of 
eig^'nvalues exceeded. Using this information, an iterative 
scheme in VIPASA identifies two values of F^^CTOR that bracket 
the desired eigenvalue. The difference between these two values 
of FACTOR can be made arbitrarily small, depending upon PASCO 
convergence criteria input CONVl and C0NV2. 

In this report, a quantity is introduced that has essentially 
the; same meaning as FACTOR. That quantity is denoted F. The 
quantities FACTOR and F differ in that whereas FACTOR is always 
the solution of an eigenvalue analysis in VIPr\,>.-. r’ is identified 
with the word FACTOR in the VIPASA printout, F may not be the 
solution of a VI*’ASA eigenvalue analysis if an adjusted shear 
analysis is used in PASCO. Otherwise, FACTOR and F are identical. 

For all analyses in PASCO, the scaler F is defined by 
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in which N , N , and N are inplane loads, P is the lateral 
X y xy 

pressure, is a bending moment, VT is a change in temperature, 

and f is a frequency. Whereas VIPASA C 2 ui have a fixed load 
system that is added to the left side of equation (39) , PASCO 
does not allow fixed loads. The entire input vector to I 

in equation (39) is scaled up or down with the quantity F to 
obtain the vector that provides the desired eigenvalue. The 
change in temperature represented by VT in equation (39) is ply 
dependent and can, therefore, be made to vary throughout the 
structure. In equation (39) , the product of F and VT 

indicates the scaling of that distribution. 

During the eigenvalue analysis, eigenvalues can be calculated 

for many values of half-wavelengths X, auid more than one 

eigenvalue can be calculated at a given value of A. There is a 

value of F associated with each of these eigenvalues. 

As an exaunple, assume that the only two nonzero eluents in 

the input vector on the left side of equation (39) are and f. 

The response of a stiffened panel might be similar to that sho%m 

in figure 20. The solid curve indicates combinations of 

and f that give the lowest eigenvalue. The value of 

that causes buckling is N ; the natural frequency of the 

*cr 

unloaded panel is f . Let the input values of N and f be < 

represented by the solid circular symbol. The dashed line that 

passes through both the origin and the circular symbol indicates | 

y 

i 

the locus of values of N and f that are considered by I 

VIPASA as possible solutions to the eigenvalue problem. The | 
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Figure 20.- Response of hypothetical stiffened panel showing 
combinations of Nx and f that provide the lowest 
eigenvalue, and a geometric interpretation of F. 

direction cosines of this line are defined by the input values 

represented by the circular symbol. The unknown is the di ;ance 

from the origin to the point at which the dashed line intersects 

the solid curve. This unknown is denoted F and can be thought 

of as the ratio of the distance AB from the origin to the solid 

curve to the distance AC from the origin to the circular symbol. 

In this hypothetical exeunple, F is approximately 0.75. Other 

examples, involving combined loads with or without vibration 

frequency, are treated in the seune manner. 

Smeared Orthotropic Stiffnesses 

Six smeared orthotropic stiffnesses are calculated by PASCO 
and are i.. luded in the output. These stiffnesses are denoted 

• Aj^j^ longitudinal extensional stiffness 

• A 22 transverse extensional stiffness 

• A ^3 shear stiffness 

• longitudinal bending stiffness 
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and D 22 calculated within 


• D 22 transverse bending stiffness 

• Djj effective twisting stiffness 

The stiffnesses ^ 22 ' ^33' 

VIPASA using the VIPASA stiffness matrix, which relates forces 
amd moments along the edges of a repeating element to the 
corresponding displacements and rotations. The stiffness matrix 
is evaluated at F * 0, X = FSTIFF • L (where FSTIFF is input, 
default =10) to approach the result obtained for a uniform 
edge loading. These stiffnesses are equivalent to the 
corresponding stiffnesses in the leuninate force - distortion 
relationships given in equations (1) and (35) . 

The stiffnesses ^11' *^33 calculated 

within VIPASA, but are, instead, calculated with formulas as 
follows : 

The smeared extensional stiffness At 
type stiffness given by 

(A ) ^^ 

'^ 12 . ' 

lA, 


11 


is defined as an ET- 






‘ 22 : 


^i 


(40) 


in which the subscript i refers to plate element i, the A^j^ 

are laminate stiffnesses defined by equation (1) , the summation 

extends over all elements in one period of the stiffened panel, 

and b„ is the width of one period, 
s 

The smeared orthotropic bending stiffnesses ^ 22 ' 

0^3 are appropriate to use in the following differential 
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equation for lateral deflection of an orthotropic plate with 
lateral loading q. 


“ii ^ 
3x^ 


9^w . pj _ 

33 . 2. 2 °22 4 ^ 

3x 3y 3y 


( 41 ) 


The smeared bending stiffness is an El-type stiffness 

given by 



in which is the distance from the centroid of the cross 

section to the centroid of plate element i, and 0 is the angle 
plate element i makes with the horizontal. 

The formula for calculating the effective twisting stiffness 
0^2 depends upon whether the panel is an open-section panel, 
such as a blade-stiffened panel, or a closed-section panel, such 
as a hat-stiffened panel. For open-section panels is 

given by 
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(43) 


°33 bg S °33^^ 


in which the summation extends over all elements in one period of 


the stiffened panel, and the D 
defined in equation (35) . The 


jk 


'12 


are laminate stiffnesses 
term is included in 


equation (43) to make correct for equation (41) . For 

f losed-section panels, is given by 


D 
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(44) 


in which A is the area enclosed by the closed section in one 
period a d the summation extends only over those elements making 
up the closed section. 


Adjusted Analysis for Shear Buckling 
Rationale for adjusted analysis approach .- With the VIP ASA 
analysis, the boundary conditions on the side edges (the edges 
parallel to the stiffeners in figure 21) can be specified and 
modeled correctly. However, the boundary conditions on the end 
edges (the edges normal to the stiffeners in figure 21) cannot 
be specified. The boundary conditions on the ends arise from the 
displacement shape assumed in equation (38) . In the case of 
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PANEL BOUNDARIES 

■ ■ ■ STIFFENERS 

NODE LINES 



Figure 21.- Node lines determined by VIPASA in the case of a 

shear loading. 

of loadings involving shear, the displacement shape and resultant 
nodal pattern produce boundary conditions on the ends that are 
not compatible with a finite rectangular panel. This incompat- 
ibility causes VIPASA to underestimate the X = L buckling 
load when the loading involves shear. 

The adjusted shear analysis is an attempt to "rectangularize" 
the nodal pattern for the A = L buckling load and, thereby, 
provide a more accurate A = L buckling analysis for loadings 
involving shear. It is assumed that node lines would be more 
compatible with a finite rectangular panel if boundary conditions 
were modeled correctly on edges normal to the stiffeners than if 
boundary conditions were modeled correctly on edges parallel to 
the stiffeners. This assumption follows from the belief that, 
in such an analysis, the stiffeners would tend to produce node 
lines that are generally parallel to the stiffeners, as shown in 
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Figure 22.- Node lines for hypothetical VIPASA shear buckling 
solution for stiffened panel rotated 90©. 

figure 22. It might appear that such an analysis could be 

carried out with VIPASA if the stiffened panel were first simply 

rotated 90°. However, VIPASA cannot solve this analysis problem 

because the VIPASA solution (eq. (38)) requires that the 

trigonometric solution be in the direction in which the panel is 

uniform, which, for a stiffened petnel, is in the stiffener 

direction. Let the value of F for this target problem be 

denoted F^ where F is defined in equation (39) , d refers 

to discrete stiffeners, and 90 refers to the panel rotated by 

90°. 

Calculation of adjusted buckling load .- Although VIPASA 
cannot calculate F^ VIPASA can solve a similar, simplified 

problem. If the stiffened panel is replaced by an equivalent 
orthotropic panel with smeared stiffnesses, the resulting panel 
is uniform in both directions. For this case, the panel can be 
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rotated 90o and boundary conditions can be modeled correctly on 
the edges normal to the stiffeners. The result of such an 
analysis would be similar to that shown in figure 22. Let the 
value of F for this smeared orthotropic panel be denoted 

^s,90* 

In the adjusted analysis approach, it is assumed that 
can be approximated by 


^d,90 



^s,90 


(45) 


where all three values of F on the right side of the equation 
are calculated with VIPASA, and the models used in the analyses 
are illustrated in figure 23. 

In equation (45) , both smeared solutions are based on the 
orthotropic stiffnesses discussed in the section entitled Smeared 
Orthotropic Stiffnesses. All other stiffnesses are assumed to be 
zero. The quantities F^ q and F^ q are calculated for X = L. 

The quantity F^ is associated with the lowest of the buckling 
loads calculated for X = W, W/2, W/3. . . W/MINLAM, where W is ♦ 

the panel width. Multiplying F^ by the ratio of 
F f. is an attempt to remove analysis inadequacies caused by 
representing the discretely stiffened panel by a smeared ortho- 
tropic panel. Note that F^ q is the standard VIP ASA solution. 

The input parameter SHEAR is used to indicate whether the 

adjusted analysis is to be used for the X = L buckling load. ; 

i 
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Figure 23.- Analysis models used to obtain adjusted 
solution for shear buckling. 

If SHEAR = 0, the standard VIPASA analysis is used. If SHEAR > 0, 
the adjusted analysis is used. When SHEAR > 0, the value of the 
twisting stiffness used in calculating the smeared orthotropic 
plate buckling load is the product of SHEAR and the value of the 
twisting stiffness calculated by equations (43) and (44). A 
value of SHEAR less than 1 is generally appropriate for a panel 
composed of closed section stiffeners, such as a hat-stiffened 
panel. 

If the adjusted analysis is selected for the X - L buckling 
load, PASCO automatically carries out the three analyses on the 
right side of equation (45) , and chooses for the adjusted 
solution the smaller of the following values of F: (1) ^*(5 go 
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calculated from equation (45) , and (2) *"0 99 calculated 

directly by VIPASA. 

To summarize the various possibilities for F: 

* When the adjusted shear analysis is used (SHEAR ^ 0 
and A = L), F is the smaller of F^ and F^ 

• For all other cases (SHEAR =0 or A L) , F is Fji q* 

The appropriate value of F is then used to calculate bending 
loads and constraints on buckling or vibration. 

The adjusted analysis is an engineering approximation, and 
engineering judgment should be used in its application. For 
example, the smeared stiffening approach must be compatible with 
the F « and F buckle mode shapes. In both cases, the 

buckle length transverse to the stiffening must be greater than 
2.5 times the stiffener spacing. If the adjusted analysis is 
used and if it is appreciably greater than the standard VIPASA 
analysis, then a factor of safety of 10 percent to 20 percent 
is recommended for the A = L buckling load. For sizing 
purposes, this factor of safety can be introduced with CLAM(l) = 
1.1 or 1.2. (See ref. 9 for additional discussion and examples.) 

Example . - An example which illustrates the approach used in 
the adjusted shear analysis is presented in figure 24. In this 
figure, buckling interaction curves for shear and compression 
are shown for a 76.2 cm (30 in.) square, blade-stiffened panel 
having six stiffeners. The desired boundary conditions are 
simple support on all four edges, a condition that cannot be met 
with VIPASA if shear is present. The four curves represent the 
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Figure 24.- Comparison of predicted buckling loads from 
various analysis models for blade-stiffened panel 
subjected to combined longitudinal compression and 
shear loadings. 

four solution approaches just discussed and are identified in the 
figure key. In particular, the solid curve represents the 
standard VIPASA analysis, and the highest curve represents the 
solution obtained using equation (45) . The circular symbols 
indicate results obtained with the STAGS computer program (ref. 10) . 
In the STAGS analysis, the panel was modeled in detail with 
discrete stiffeners, and the desired simple support boundary 
conditions were maintained on all four edges. For this problem, 
the standard VIPASA analysis greatly underestimates the shear 
buckling load. Either of the two upper curves provides a 
reasonably accurate estimate of the correct result obtained with 
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the general two-dimensional STAGS analysis. As explained earlier, 
if an adjusted analysis were desired, PASCO would automatically 
choose the lower of the two upper curves. (See examples, ref. 9.) 


SIZING 

The computerized structural sizing approach used in PASCO 
is based on nonlinear mathematical programming techniques. Sizing 
variables are automatically adjusted to obtain a design that 
minimizes an objective function subject to a set of inequality 
constraints. Approximate analysis techniques are used to improve 
computational efficiency. 
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Problem Statement 

The general problem statement is: find values for the set 

of variables to 

• Minimize an obj -tive function OBJ (X^^) 

• Subject to 

• Behavioral constraints: G^ (X^) £ 0 

• Side constraints: VLB^^ £ VUB^ 

where 

X^ are the sizing variables 

VLB^ are the lower bounds on the sizing variables 
VUB^ are the upper bounds on the sizing variables 
VIPASA and other analyses are used to evaluate the constraints 
Gj. CONMIN (refs. 11 and 12) is used to solve the resulting 
mathematical progreunming problem. 
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Sizing Variables 


In PASCO, the sizing variables are the plate element widths, 
denoted b, the ply thicknesses, denoted t, and the ply 
orientation angles, denoted 6. Any set of widths, thicknesses, 
and orientation angles can be selected as the active sizing 
variables. The remaining widths, thicknesses, and orientation 
angles can be held fixed or linked linearly to the active sizing 
variables. Upper and lower bounds can be specified for the 
sizing variables. 


Objective Function 

W/A 

The objective function is the panel mass index — ^ , 
the panel mass per unit area divided by the panel length. This 
is the quantity denoted OBJ in CONMIIn. The area A is the panel 
planform area shown in figure 25. Since the panel length L is 
fixed, the quantity that is minimized becomes the panel mass per 
unit width. 



m 

A - LENGTH- WIDTH 

Figure 25.- Panel planform area A. 
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Constraints 


Constraints are inequality requirements that must be met 
during sizing to provide an acceptable design. In addition to 
upper and lower bounds on the sizing variables, denoted side 
constraints, there are behavioral constraints on buckling, 
material strength, stiffness, and vibration frequency. CONMIN 
requires that these constraints be written in the form 


G. 


3 


< 0 


(46) 


In PASCO, the constraints =re normalized in order that all 
constraints be of the same order of magnitude. The specific 
forms for the constraints are given xn the following sections. 

Buckling or vibration .- Constraints on the buckling load 
or vibration frequency can be written in the form 


F(X.) 

G = 1 J 

CLAM(Xj) 


(47) 


in which F is defined in equation (39) , and CLAM can be used to 
specify a margin of safety at specific wavelengths. There can 
be simultaneous buckling or frequency constraints for many values 
of X, and there can be many buckling or frequency constraints 
for each value of X . 

In the coding within PASCO, F is replaced by N /N , 

cr *input 

which is equivalent to F. If the adjusted shear analysis is 
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selected, the appropriate analysis is used to compute N . If 

*cr 

N is zero, a small positive value of N is 

input input 

automatically introduced within PASCO. 

Material strength .- Three material strength criteria are 

built into PASCO: maximum lamina stress, maximum Icunina 

mechanical strain, and Tsai-Wu (ref. 13) . In the maximum stress 

criterion, tension and compression limits are placed on 

and in each lamir . The maximum lamina mechanical strain 

criterion is defined similarly, except that the thermal strain is 

subtracted from the total strain to provide the mechanical strain 

in each element . For the maximum lamina stress and maximum 

lamina mechanical strain criteria, the material strength 

constraint is written in the form 


G = ^ 1 M8) 

allow 

in which S is a lamina stress or mechanical strain, and S ,, 

' allow 

is the corresponding maximum allowable value. The input quantity 
ALLOW is used to prescribe the allowable values used in the 
» material strength criteria. 

In the Tsai-Wu criterion, the stress state is defined by 


4 ) - + F22O2 


^ 66'^12 ^^ 12 '^ 1°2 


(49) 
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where Pj^, Pj# Pj^j^» P 22 / Pgg * 1 ^® automatically calcu* 

lated from the allowable stresses that are included in the input 
ALLOW, and P ^2 included in the input ALLOW. The Tsai-Wu 
strength constraint is defined as 


G = (^ - 1 


(50) 


The user may incorporate his own material strength criteria 
by writing additional subroutines. 

Stiffness.- Stiffness constraints are written as 


G = 1 


Stiffness 

Stiffness lower limit 


(51) 


and 


Stiffness 

Stiffness upper limit 


(52) 


in which the stiffnesses that can be constrained are: 

• extensional stiffness 

• A ^2 shear stiffness 

• bending stiffness 

These stiffnesses are "smeared" orthotropic stiffnesses for Lhe 
overall panel, not individual plate element stiffnesses. 
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Approximate Analysis 

The approximate analysis approach used in PASCO is depicted 
in figure 26. It is similar to the approach proposed in 
reference 14. The procedure consists, conceptually, of three 
modules: an analysis module, a Taylor series module, and a 

sizing module. 


ANALYSIS 

MODULE 


TAYLOR SERIES 
MODULE 


RESIZING 

MODULE 



Figure 26.- General approach used in PASCO. 

Analysis module .- In the analysis module, all constraints 
are calculated with VIPASA and supporting subroutines . The program 
identifies the critical constraints and, using a two-point forward 
difference approximation, calculates the derivatives of the 
critical constraints with respect to the sizing variables. The 
values of the constraints and derivatives are then passed to the 
second module, the Taylor series module. The techniques used to 
■ ?ntify critical ccnctraints are discussed subsequently. 

Taylor series module .- The Taylor series module generates a 
first order Taylor series expansion of each constraint. Expan- 
sions are of the form 
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( 53 ) 


G(X.) = G(X^) - X^) 

1 ' - Xj 

in which are the sizing variables and are the values of 

the sizing variables at the initial point of the expansion. The 
Taylor series approximations provide a reasonably accurate and 
simple representation of the constraints in the neighborhood of 
the initial point of the expansion. The Taylor series expansions 
are updated periodically to insure their adequacy. The second 
module also evaluates the objective function. 

Sizing module .- The third module contains the optimizer 
CONMIN. During sizing, the optimizer interacts only with the 
second module which contains approximate, explicit functions 
for the coiistraints and a simple expression for the objective 
function. Such as approach greatly improves computational 
efficiency. 

Sizing strategy .- The overall sizing strategy is depicted in 
more detail in figure 27. The strategy consists of a series of 
sizing cycles in which the optimizer adjusts the values of the 
sizing variables based on approximate values of the constraints 
(eq. (53)). An upper limit is imposed on the change of each 
sizing variable during each sizing cycle to insure the adequacy 
of both the list of constraints that are considered to be critical 
and the Taylor series expansions of those constraints. These 
limits to the changes in the sizing variables, referred to as 
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COMPLETE ANALYSIS 


MOVE LIMITS FOR 
n'f* SIZING CYCLE 


FIRST SIZING VARIABLE. Xi 


Figure 27.- Sizing strategy for approximate analysis, shown 
in two-sizing-variable space. 

move limits, are governed by input and are indicated by the 
dashed rectangles in figure 27. (Move limits are discussed 
in the next section.) The solid circular symbol at the center of 
each rectangle in figure 27 represents the point at which the 
Taylor series expansions are carried out for each sizing cycle. 

The end point of one sizing cycle becomes the initial point of the 
next sizing cycle. Accurate values of the constraints and deriva- 
tives of the constraints are then recalculated, and new Taylor 
series expansions are generated. Ten sizing cycles are usually 
adequate to obtain convergence if the initial design is reasonably 
well chosen. The number of sizing cycles is controlled by the 
input parameter MAXJJJ and not by any convergence criterion. 
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Move limits .- The move limits that are generated internally 


for each sizing cycle are given by 


VLBj = if. 

- DVMOV^ * (SFACTR)”"^ 

• X 

i,init 

(54) 


VUBj . iT 

+ DVMOV^ • (SFACTR)*'”^ 

* X 

i,init 

(55) 



where 

• VLB^ and VUB^ are the sizing variable lower and upper 
bounds used by CONMIN in a sizing cycle 

• are the values of the sizing variables at the 
beginning of a sizing cycle 

• DVMOV^ is an input vector 

• SFACTR is an input scaler 

• n is the sizing cycle number 

• X. are the initial (input) values of X. 

One of the objectives of equations (54) and (55) is to reduce 
the move limits as the sizing progresses. Overall lower and 
upper bounds on the sizing variables override the lower and 
upper boiinds for a sizing cycle calculated in equations (54) 
and (55). Values of DVMOV =0.2 and SFACTR =0.8 generally 
provide reasonable answers. 

Identifying critical buckling and frequency constraints . - 
For simplicity, buclcling terminology rather than eigenvalue 
terminology is used to describe the logic for identifying 
critical eigenvalue constraints for the Taylor series module. 
However, the discussion also applies to the frequency constraints. 
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The critical buckling constraints are identified by buckling 
h?lf-wavelength A. Selecting these critical values of X is 
a multistep process which begins by constructing a table of 
potentially critical values of X. This table always contains 
X = L. The tcible also contains values of X specified in the 
input NLAM. Also added to the table is each value of X for 
which the buckling load meets both of the following two require- 
ments . 

• The buckling load is a relative minimum (X = L/8 in 
figure 17) , and 

• The buckling load is lower than the buckling load for the 
preceding value of X in the X table. The X table 
is ordered according to decreasing values of X . 

As the sizing progresses, new values of X meeting these two 
requirements are identified and added to the table. 

Buckling constraints are retained for a maximum of MAXL 
values of X, in which MAXL is an input paraimeter. If there are 
more than MAXL values of X in the X table, logic is included 
which divides the range of m values (m = 1, 2, 3,..., MINLAM) 
into regions and retains the most critical constraint (s) in each 
region. The larger MAXL, the larger the number of regions. At 
this stage, X = L is always retained in the table of potentially 
critical values of X. 
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The number of buckling constraints can be further reduced 
with the parameter GRANGE. All buckling constraints for which 
F/CLAM > GRANGE are eliminated. For example, for GRANGE = 2 
any buckling constraint for which G is less than -1.0 is 
eliminated. To insure computational efficiency, all buckling 
constraint elimination is based on the rough calculation of the 
buckling loads where the convergence criterion is CONVl. Fine 
calculations of the buckling loads and calculation of the 
derivatives of the buckling loads are not carried out for buckling 
constraints that are eliminated. 

Identifying other critical constraints .- For constraints 
other than buckling, GRANGE is the sole constraint deletion 
mechanism. For constraints involving a lower bound (buckling, 
stiffness, frequency) the minimum value of G retained is 


G = 1 - GRANGE 


(56) 


For cc .straints involving an upper bound, (stiffness, material 
strength) the minimum value of G retained is 


G 


1 

GRANGE 


1 


(57) 


Calculation of Derivatives of Buckling Loads 
The derivatives of the buckling loads with respect to sizing 
variable are calculated using the following numerical 

approximation 
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(58) 


dF 

dX. 

1 


» 


F(X^ + AX^) 
AX- 

i 


- F(X^) 


Derivatives are calculated only for those values of the half- 

* wavelength \ identified as being critical. Two methods are 

available in PASCO for calculating the perturbed solution 

F(X^ + AX^) . One method uses the same general approach as that 

used to calculate the nominal solution F(X.). The other method 

1 

uses a much faster approximate technique. The user can select 
the method with the input parameter JDER. 

In the first method, F(X^ + AX^) is calculated using an 
iterative technique that is the same as that used to obtain the 
nominal solution F(X^). The number of iterations required to 
obtain the perturbed solution is, however, reduced somewhat by 
restricting the solution for F(X^ + AX^) to a narrow band 
centered* on the nominal solution F(X^) as shown in figure 28. 

The solid curve indicates the value of the buckling determinant 

as a function of F for the nominal case. The dashed curve 

I 

? gives the same information for the perturbed case. The nominal 

f 

t 

^ solution, the band width, and the perturbed solution are 

indicated in the figure. 

* In the second method, the perturbed solution F(X^ + AX^^) is 

estimated using an approximate technique illustrated in 
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NOMINAL CASE 


BUCKLI:<IG 

DETERMINANT 



I 


Figure 28.- Buckling determinant as a function of F for 

nominal values of the sizing variables and for a perturbed 
sizing variable. 

3 

figure 29. The approximate method consists of the following 
steps: 

• Using the value of F obtained from the nominal solution 
shown at point 1, the value of the determinant is 
calculated for the perturbed case shown at point 2. 

• The slope of the dashed curve at point 2 is assumed to be 
the same as the slope of the solid curve at point 1. 

• The perturbed solution at point 3 is estimated using a 
linear approximation from point 2. 


3 

The authors are indebted to Prof. Fred W. Williams, 
University of Wales, Institute of Science and Technology, for 
suggesting this technique. 
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Figure 29.- Approximate method for calculating perturbed 

buckling load. 

In all cases, the iteration scheme used to converge on an eigen- 
value at point 1 involves an interval halving strategy. If 
certain numerical conditions are met during iteration, a linear 
interpolation strategy is introduced. Consider the two values 
of F that bracket the desired eigenvalue when the linear 
interpolation strategy is introduced. If these two values of F 
differ by at least five percent, then the second method can be 
used to calculate F(X^^ + Ax^) . In this case, the slope at 
point 1 is taken to be the last slope calculated in the linear 
interpolation strategy. If the numerical characteristics of the 
problem are such that interpolation is not used, or if, when 
interpolation is introduced, the two values of F that bracket 
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the solution differ by less than five percent, the first method is 
automatically used to calculate P(X^ + AX^) . 

In the case of the adjusted analysis for shear buckling, the 
derivative of the adjusted buckling load for A * r, is taken to be 


dF _ _F _±_ „ 

dX, " F , - dX. ^d,0 
1 a, u 1 


(59) 


where the derivative of F, f. 

a, u 

approximation in equation (58) 
^s,90* 


is calculated using the numerical 
, and F is the smaller of F^ 


Multiple Load Conditions 

PASCO can treat sizing problems with multiple loading 
conditions. This means that panels are sized to meet the design 
requirements (constraints) for load condition 1 and load condi- 
tion 2 and load condition 3, etc. The number of allowable load 
conditions is large. For limitations, see reference 1. 

Using PASCO notation, quantities that can depend upon the 
load condition are the inplane design loads NX, NY, and NXY; 

( 

the lateral pressure PRESS; the bow ECC; the change in temperature 
TEM; the bending moment MX; the design frequency FREQ; the 
material properties El, E2, E12, ANUl, RHO, ALFAl, and 

ALFA2; the stiffness requirements AllL, AllU, A33L, A33U, 

DHL, and DllU; and the allowables ALLON used for the material 
strength criteria. 
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sizing Example 


An example which focuses on the buckling constr . t aspects 
of the sizing code is presented in figures 30 through 3"?. In 
this example, the panel whose buckling response was shown in 
figure 17 is sized to meet a buckling requirement of = 700 kN/m 
(4000 Ibf/in) . No other constraints or loadings are imposed. The 
panel is assumed to be perfect. 



L 


Figure 30.- Longitudinal buckling load as a function of 

buckling half-wavelength for two blade-stiffened panels: 
the initial design and the final design. 
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Figure 31.- Buckling loads as a function of sizing cycle number. 

The buckling response for the initial design and the buckling 
response for the final design are both shown in figure 30. As 
can be seen in the figure / the initial design buckles at about 
70 percent of the design load at buckling half-wa' elengths of 
X = L and X = L/8. The final design meets the buckling 
requirement for all values of X. 

The buckling load history for this sizing exeunplc is shown 
in figure 31. For the initial design, the critical buckling half- 
wavelength identified by the analysis module is X = L. The half- 
wavelength X - L/8 is not identified because the buckling load 
for X = l/ 8 is slightly higher than the buckling lo^.d for 
X = L. As the sizing progresses, other values of X are 
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Figure 32.- Panel mass index as a function of sizing cycle number, 
identified and ar^' added to the A table. For example, at the 
second sizing cycle, X = L/7 is identified, and at the third 
sizing cycle L/8 is identified. 

The manner in which the mass index varies during the sizing 
is shown in figure 32. During the first cycle, the only buckling 
constraint available to PASCO was for the A = L mode- The code 
increased that buckling load while decreasing the panel mass. At 
the beginning of the second sizing cycle, PASCO identified the 
A = L/7 mode and, during that cycle, added materi'i to the panel 
to increase that buckling load. the beginning of the fourth 

sizing cycle, all buckling constraints were satisfied. During 
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fourth and subsequent sizing cycles, the mass was reduced while 
the buckling strength of the panel was maintained. 

In the above example, the first sizing cycle was counter- 
productive because the only buckling constraint available to 
PASCO was for the X = L mode. Convergence would have been 
improved if a local buckling constraint had been made available 
with the input NLAM. In this case, a good choice would have been 
NLAM = 7 or 8. 

CONCLUDING REMARKS 

This report has discussed certain aspects of the computer 
code denoted PASCO, which can be used for analyzing and sizing 
uniaxially-stiffened ccmposite structural panels having a general 
configuration . 

In PASCO, buckling loads, lamina stresses and strains, smeared 
orthotropic stiffnesses, and vibration frequencies can be calcu- 
lated for a variety of typical loading conditions. These same 
quantities can also be used as design requirements during sizing. 
Sizing is based on nonlinear mathematical programming techniques 
ir. which the mass of the panel is minimized subject to satisfaction 
of the design requirements. 
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